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I. INTRODUCTION 



Quantum systems are characterized by a number of discrete states the structure of which 
is more or less comphcated. These quantum systems do, however, not exist isolated from 
other systems. Most of them are embedded in an environment, e.g. in the continuum of 
decay channels. System and environment interact with one another, and it is this interaction 
that allows us to study the properties of the system. The feedback of the interaction with 
the environment onto the properties of the system itself is an old problem raised in the very 
beginning of quantum mechanics. It becomes the more important the smaller the system is. 

As an example, most states of a nucleus are embedded in the contiuum of decay channels 
due to which they get a finite life time. In other words: the discrete states of a nucleus shade 
off into resonance states with complex energies £k = Ek — ^Tk- The Ek give the positions in 
energy of the resonance states while the widths are characteristic of their life times. The 
Ek may be different from the energies of the discrete states, and the widths may be large 
corresponding to a short hfe time. Nevertheless, there is a well defined relation between 
the discrete states characterizing the closed system, and the resonance states appearing in 
the open system. The main difference in the theoretical description of quantum systems 
without and with coupling to an environment is that the function space of the system is 
supposed to be complete in the first case while this is not so in the second case. Accordingly, 
the Hamilton operator is Hermitian in the first case, and the eigenvalues are discrete. The 
resonance states, however, characterize a subsystem described by a non- Hermitian Hamilton 
operator with complex eigenvalues. The function space containing everything consists, in 
the second case, of system plus environment. 

The mathematical formulation of this problem goes back to Feshbach [1] who introduced 
the two subspaces Q and P, with Q + P — 1, containing the discrete and scattering states, 
respectively. Feshbach was able to formulate a unified description of nuclear reactions with 
direct processes in the short-time scale and compound nucleus processes in the long-time 
scale. Due to the high excitation energy and high level density in compound nuclei, he intro- 
duced statistical approximations in order to describe the discrete states of the Q subspace. 
A unified description of nuclear structure and nuclear reaction aspects is much more compli- 
cated and became possible only at the end of the last century (see [2] for a recent review) . In 
this formulation, the states of both subspaces are described with the same accuracy. All the 
coupling matrix elements between different discrete states, different scattering states as well 
as between discrete and scattering states have to be calculated in order to get results that 
can be compared with experimental data. This method has been applied to the description 
of light nuclei by using the shell model approach for the discrete many-particle states of the 
Q subspace [2]. 

In the unified description of structure and reaction aspects, the system is described by an 
effective Hamiltonian that consists of two terms: the Hamiltonian matrix H of the closed 
system with discrete eigenstates and the coupling matrix between system and environment. 
The last term is responsible for the finite lifetime of the resonance states. The eigenvalues 
of Ti. are complex and give the poles of the S matrix. The motion of these eigenvalues as 
a function of a certain parameter is discussed in many papers (see the recent review [2]). 
The statistics of complex eigenvalues and the corresponding nonorthogonal eigenvectors for 
non-Hermitian random matrices are recently considered in [3]. 
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The coupling matrix elements W^' between the resonance states and the continuum are 
seldom studied. Their relation to the coupling matrix elements Wj:'^ between the corre- 
sponding discrete states and the continuum can be expressed by the mixing coefficients that 
appear in the representation of the wave functions of the resonance states in the set of wave 
functions of the discrete states. Generically, the relation between the wave functions of the 
resonance states and those of the discrete states is complicated since the number of discrete 
states of a realistic system is large. Many of them can contribute to the wave function of a 
certain resonance state, almost independently of their energetical distance, see e.g. [4]. The 
numerical results for the coupling matrix elements W^*^ of nuclear states show a non-trivial 
energy dependence, especially at high level density [2] . In the statistical approach to nuclear 
reactions and the application of this approach to some other reactions they are, however, 
assumed to be simple, energy-independent parameters such as the Wj:'^ , e.g. [5-9]. 

It is the aim of the present paper (i) to study the energy dependence of the coupling 
matrix elements W^'^ between resonance states and continuum that follows immediately 
from the unitarity of the S matrix, and (ii) to study the behaviour of the wave functions of 
the resonance states in the overlapping regime since they determine the energy dependence 
of the Wj:'^ in numerical calculations performed for special systems. Most interesting is the 
behaviour of the wave functions near avoided level crossings in the complex plane. 

As the results show, the coupling coefficients Wj:" are, generically, energy dependent. The 
energy dependence is, however, not important as long as the distance in energy between the 
resonance states is larger than the sum of their widths. This result is in full agreement 
with the statement of the authors of the review [5] who restricted the application of their 
approach to the non- overlapping regime. In the overlapping regime, however, the energy 
dependence of the coupling coefficients Wj:'^ can not be neglected. It is a consequence of 
the unitarity of the S matrix and causes nonlinear terms in the S matrix at high level 
density. Nevertheless, the fine shape of the resonances can equivalently be described by the 
energy independent Wj^'^ in many cases. The Wj^'^ lose, however, their physical meaning 
in the overlapping regime [10]. As to the wave functions of two resonance states k and / 
in the neighbourhood of an avoided level crossing in the complex plane, they are mixed: 
Pk^k =1= ^ A^i- The corresponding phase changes in approaching the critical value of the 
parameter at which the levels avoid crossing (where (3k — ^Pi), are caused by nonlinear terms 
in the Schrodinger equation. These terms are, finally, responsible for the energy dependence 
of the coupling coefficients W^'^ between system and continuum in numerical calculations. 
This result following from the behaviour of the wave functions of the resonance states at 
high level density, coincides with that following from the unitarity of the S matrix. 

The paper is organized in the following manner. In Sect. II, the main ingredients of 
the unified description of structure and reaction aspects of a quantum system embedded 
in a continuum are given. These are the effective Hamiltonian Ti. of the system and the 5" 
matrix which both are derived by solving the Schrodinger equation in the full function space 
with discrete and continuous states. Further, some properties of the spectroscopic values 
that characterize the system are sketched with special emphasis of their behaviour in the 
overlapping regime. In Sect. Ill, the coupling coefficients W between system and continuum 
are directly obtained in the one-channel case by starting from a unitary S matrix. The 
nonlinear effects appearing in the overlapping regime are discussed. The wave functions 
at avoided level crossings in the complex plane and in its neighbourhood are derived in 
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Sect. IV. They are energy dependent and their phases change in a certain range aroTind the 
critical value of the parameter at which the resonance states avoid crossing. Sect. V contains 
concluding remarks on the energy dependence of the coupling coefficients W^'^ appearing in 
a model with unified description of structure and reaction aspects at high level density. 



II. EFFECTIVE HAMILTONIAN AND S MATRIX FOR A QUANTUM SYSTEM 

EMBEDDED IN A CONTINUUM 



In the unified description of structure and reaction aspects of quantum systems, the 
Schrodinger equation 

(i/fuii _ E)xif{E) = (1) 

is solved in a function space containing everything, i.e. discrete as well as continuous states. 
The Hamilton operator H^'^^^ is Hermitian, the wave functions ^ depend on energy as well 

as on the decay channels and all the resonance states of the system. Knowing the wave 
functions "^{E), an expression for the S matrix can be derived that holds true also in the 
overlapping regime, see the recent review [2]. It reads 



gi(<5c-5^/) 



iT q(1) q(2) 

Occ' iJcc' '-'cc' 



(2) 



where 5c is the phase shift in channel c, S),J is the smooth direct reaction part related to 
the short-time scale, and 

^(2)^.^ lili (3) 

k=l — J^k + 2^ k 

is the resonance reaction part related to the long-time scale. Here, the Sk — — ^ ^-re 
the complex eigenvalues of the non-Hermitian Hamilton operator 

= Hqq = Hqq + Hqp HpQ (4) 

appearing effectively in the system {Q subspace) after embedding it into the continuum (P 
subspace). It is iJ^"" = Hqq + Hqp+Hpq + Hpp where Hqq is the Hamiltonian of the closed 
system and Hpp is that for the environment (scattering states). The two terms HpQ and 
Hqp characterize the coupling between the two subspaces. They appear in the source terms 
of the equations in either subspace when Q + P — 1 as well as in the effective Hamiltonians 
of the subspaces. The effective Hamiltonian in the Q subspace is given in (4), and an 
analogous expression for the effective Hamiltonian in the P subspace can be written down, 
see [2]. Usually, Re (HqpG^^^ Hpq) ^ and Re {H) = Hqq + Re {HqpG^^^ Hpq) ^ Hqq, 
see [2] for a detailed discussion. The G^'' in (4) are the Green functions in the P subspace, 
and Im {Hqp G^p^ Hpq) is determined by the coupling matrix elements 



7fe = {CE\V\^k) (5) 

between the resonance states and the scattering states. The £k are energy dependent func- 
tions. They arc complex and also energy dependent functions. In (5), the are the scatter- 
ing wave functions and the are the eigenfunctions of 7i, Eq. (4). The wave functions flk 
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of the resonance states are related to the eigenfunctions of by a Lippmann-Schwinger 
hke relation [2] , 

Qk = {l + G^^^HpQ)^k. (6) 
The eigenfunctions of 7i are bi-orthogonal, 

{^*\^,)^S,i (7) 

so that 

(l>fe|l>fc) = Re((l>fe|l>fc)) ; Ak = {h\^k) > 1 (8) 

and 

{h\^i^k) = tM{^k\^i^k)) = -{^iM^k) ; Bi^'' = \{h\^i^k)\>o . (9) 

It should be noticed that the standard normalization = Ski is equivalently to Eqs. 

(7) up to (9) for all k, I but those with Sk = Si (double pole of the 5* matrix) where — > oo 
and BI ^ oo. As a consequence of (8) holds [2] 

f, = ^4^ < E\r,\' . (10) 

c 

The main difference to the standard theory is that the F^, 7^ and are not numbers 
but energy dependent functions [2]. The energy dependence of Im{£k} — — | Ffc is large near 
to the threshold for opening the first decay channel. This causes not only deviations from 
the Brcit Wigner line shape of isolated resonances lying near to the threshold, but also an 
interference with the above-threshold " tail" of bound states [2] . Also an inelastic threshold 
may have an influence on the line shape of a resonance when the resonance lies near to the 
threshold and is coupled strongly to the channel which opens [11]. Also in this case, Tf. 
depends strongly on energy. In the cross section, a cusp may appear in the cross section 
instead of a resonance of Breit Wigner shape. Both types of threshold effects in the line 
shape of resonances can explain experimental data known in nuclear physics [2] . They can 
not be simulated by a parameter in the S matrix. 

The energy dependence of Ek and may be important also far from decay thresholds [2] . 
Characteristic of the motion of the poles of the S matrix as a function of a certain parameter 
(which may be also the energy E of the system) are the following generic results obtained 
for very different systems in the overlapping regime: the trajectories of the S matrix poles 
avoid crossing with the only exception of exact crossing when the S matrix has a double (or 
multiple) pole. At the avoided crossing, either level repulsion or level attraction occurs. The 
flrst case is caused by a predominantly real interaction between the crossing states and is 
accompanied by the tendency to form a uniform time scale of the system. Level attraction 
occurs, however, when the interaction is dominated by its imaginary part arising from the 
coupling via the continuum. It is accompanied by the formation of different time scales in 
the system: while some of the states decouple more or less completely from the continuum 
and become long-lived (trapped), a few of the states become short-lived and wrap the long- 
lived ones in the cross section. The dynamics of quantum systems at high level density is 
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determined by the interplay of these two opposite tendencies. For a more detailed discussion 

sec [2]. 

One of these tendencies, the phenomenon of resonance trapping 

N ^ K ^ N ^ 

J2^k^ Effc ; E ffe^^o, (11) 

k=l K=l k=K+l 

appears only in the overlapping regime. It is caused by lm(7f) and means almost complete 
decoupling oi N — K resonance states from the continuum while K of them become short- 
hved [2] . Usually, K <^ N — K. The long-lived resonance states in the overlapping regime 
appear often to be well isolated from one another [10]. The few short-lived resonance states 
determine the evolution of the system. That means, quick direct reaction processes may 
appear, at large overall coupling strength, from slow resonance processes by means of the 
resonance trapping phenomenon. Meanwhile, the phenomenon of resonance trapping has 
been proven experimentally on a microwave cavity as a function of the degree of opening of 
the cavity to an attached lead [12]. In this experiment, the parameter varied is the overall 
coupling strength between discrete and scattering states. Resonance trapping may appear, 
however, as function of any parameter [2]. 

In any case, the energies and widths of the resonance states follow from the solutions of 
the fixed-point equations : 

Ek = Ek{E=Ek) and Tk = Vk{E=Ek) , (12) 

on condition that the two subspaces are defined adequately [2]. The values Ek and 
correspond to the standard spectroscopic observables. The wave functions of the resonance 
states are defined by the functions Eq. (6), at the energy E = Ek- The partial widths are 
related to the coupling matrix elements (7^)^ that are calculated independently by means 
of the eigenfunctions of H. For isolated resonances, Ak — I according to (8), and the 
standard relation Ffc = Z^dTfeT follows from (10). In the overlapping regime, the partial 
widths lose their physical meaning, since A > 1. Both functions, (7^)^ and Fj^, may even 
show a different energy dependence [2]. 

It follows immediately from (4), that the coupling of the resonance states via the con- 
tinuum induces additional correlations between the states. These correlations are described 
by the term HgpGp'^HpQ of the effective Hamiltonian H. The real part Re(HQpGp'^HpQ) 
causes level repulsion in energy and is accompanied by the tendency to form a uniform time 
scale in the system. In contrast to this behaviour, the imaginary part lm{HQpG^p^ Hpg) 
causes different time scales in the system and is accompanied by level attraction in energy. 
Thus, an essential part of the energy dependence of the eigenvalues of H is caused by the 
additional correlations of the states via the continuum. They are important, above all, at 
high level density. 

For isolated resonance states, the additional shift in energy is usually taken into ac- 
count by simulating Re(7i) = Hqq + Re {NqpOp'^ Hpg) by Hq + V, where V is assumed 
to be describable by two-body effective residual forces. It should be mentioned, how- 
ever, that Re{HQpG^p^ Hpq) can not completely be simulated by an additional contribu- 
tion to the residual two-body interaction even in the case of well isolated resonances, since 
it contains many-body effects, as follows from the analytical structure of HqpG^p^ HpQ. 
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Re{HQpGp Hpq) is an integral over energy and depends explicitly on the energies ec at 
which the channels c open. For details see [2]. 

Spectroscopic studies in the overlapping regime are more complicated. The wave func- 
tions $fc may be represented in the set of eigenfunctions of the Hermitian Hamilton 
operator H = Hqq , 

$/c = E hi . (13) 
I 

The $fc are real, while the are complex and energy dependent in the overlapping regime. 
The coefficients bki and the 7^ are complex and energy dependent, too. In this regime, the 
differences between Ti. and H can therefore not be simulated in a simple manner. Even the 
positions of the peaks in the cross section do, generally, not appear at the energies when 
the resonance states overlap [2]. They are the result of interferences between the resonance 
states. 

It should be underlined here that the expression (3) for the resonance reaction part of 
the S matrix is derived from the Schrodinger equation (1) by rewriting it in a consistent 

manner. Here, the eigenvalues £k = -E-fe — f Tfc of the effective Hamiltonian Ti, Eq. (4), as 
well as the coupling matrix elements 7^ are energy dependent functions, and the unitarity 
of the S matrix is garanteed. 

Furthermore, the different ^k{E=Ek) are neither strictly orthogonal nor bi-orthogonal 
since the bi-orthogonality relation (7) holds only when the energies of both states k and I are 
equal. The spectroscopic studies on resonance states are performed therefore with the wave 
functions being only approximately bi-orthogonal. The deviations from the bi-orthogonality 
relation (7) are, however, small as a rule. This drawback of the spectroscopic studies of 
resonance states has to be contrasted with the advantage it has for the study of observable 
values: the S matrix and therefore the cross section are calculated with the resonance wave 
functions being strictly bi-orthogonal at every energy E of the system. Furthermore, the 
full energy dependence of Ek, F^ and, above all, of the coupling matrix elements 7^ is taken 
into account in the S matrix and therefore in all calculations for observable values. 



III. UNITARITY OF THE S MATRIX 

The Breit Wigner one-level formula for nuclear reactions describes the reaction cross 
section with isolated resonances. The S matrix elements for this case read 

S^, = 1 - ^ V 'J^ — (14) 

k ^ ~ -t^k + 2 '- k 

where Ek and Ffc are the energies and widths of the resonance states k and Wk"' = 7^7^'. 
The ^k the partial widths of the states k relative to the channel c. The values Ek, Tk 
and 7^ are numbers characterizing the properties of the resonance states k. Since they are 
energy-independent values, the decay follows an exponential law. 
For an isolated resonance state k = 1 coupled to one channel, 

S^l- i ' . . (15) 
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in the energy range £"1 — ^ Fi < < £"1 + ^ Fi, and Wi = Fi due to the unitarity of the 
S matrix. The last relation follows immediately from (15) that, in the one-resonance-one- 
channel case, can be written as 

„ E — El — 7: Fi , . 

S = J ^ 16 

£ - £1 + I Fi 

when Wi = Fi. The S matrix (16) is unitary- 
Let us now consider the unitary representation of the S matrix in the one-channel case 
with two resonance states, 

{E-Ei-^Ti) {E-E^-^t^) 
{E-Ei+ ifi) {E-E2+ ffs) ■ 

Prom this expression, a possible form of the pole representation of the S matrix can be 
derived. 



£ - £1 + I Fi £ - £2 + i {E-Ei + iFi)(£ - E^ + fFa) 

1 - — ~ . ^ hri + 




2E-{Ei + E2)+ |(Fi + F2), 

TiFs \ 

2E-{Ei + E2)+ |(fi + f2); 



(18) 



It follows 



s = i -t — i~ (19) 

fe=i,2 E - Ek+ I Ffc 

in complete analogy to (14), with 

Wk = Tk(l-i = £ ^1 (20) 

H 2E-{Ek + Ei)+ ^iFk + Ti) J ^ ' 

and k,l = 1,2. I ^ k. These equations show that the coupling coefficients Wk arc complex 
and energy dependent, that Wk has a resonance behaviour at the energy {Ek + Ei)/2 with 
the width (F^ -|- Ti)/2, and that the energy dependence of the two values Wk and Ffc is 
different in the overlapping regime. In the energy region of the resonance behaviour of Wk 
caused by a neighboured resonance state I, the S matrix contains terms being nonlinear in 
energy. 

When the widths of the two states are equal, Wk — with E — > {Ek + Ei)/2. At large 
distance E Ei + E2, follows Wk — > Ffc- In this case, the two resonance states behave as 
isolated ones. When the positions Ek, Ei of the two resonance states are outside of resonance 
region of Wk and Wi, they can also be considered, to a good approximation, as isolated, and 
Wk ^ Tk, Wi ^ fi. 
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The nonlinear term creates some deviation in the resonance hne shape from the hnear 
Breit-Wigner one. This can be seen best in the case when the S matrix has a double pole, 
i.e. El = E2 = Ed and Fi = r2 = r^. In such a case, the S matrix (19) reads 

5 = l-2z. 



E-E.+ 



f 

^l-2t ^ . - — . (21) 

E-Ea+^Fa (E-Ea+^rd^ 

The second term corresponds to the usual linear term while the third term is quadratic 
(see [13]). The interference between these two parts has been illustrated in [14], where the 
cross section is shown for the case of two resonance states coupled to one channel. The 
energies and widths of the two resonance states are the same, creating a double pole of the 
S matrix. The asymmetry of the line shape of both peaks in the cross section agrees with 
(21). A similar picture has been obtained in, e.g., laser induced continuum structures in 
atoms with a double pole of the S matrix [15,16], in atom-surface collisions [17], transmission 
in quantum scattering systems [18], in a double barrier potential [19], a double-square- well 
system [20] and in a toy model for the conductance through a small quantum dot [21]. 

According to (20), the asymmetry of narrow resonances is usually larger than that of 
broad resonances: when Fi ^ it follows Wi ~ Fi while the corrections from Fi to W2 
can mostly not be neglected. In any case, the nonlinear term in Wk, eq. (20), causes a 
non-exponential decay of the two resonance states. Only when the line shape of a certain 
resonance k is of Breit-Wigner type and F^ is almost constant in a large energy region around 
Ek, the states k will decay according to an exponential law [13]. 

When the two resonance states lie at the same energy, but their widths are different, 
then follows from (19) and (20) that the contributions from both resonance states annihilate 
each other at E — Ei — E2, i.e. the cross section vanishes at that energy where the two 
resonance states lie. This destructive interference has been traced numerically in [14] for 
two resonance states coupled to one channel by varying the coupling strength between the 
states and the continuum. When Fi ^ F2, the narrow resonance appears as a dip in the 
cross section that is determined mainly by the broad resonance. This is in accordance with 
(20): for E ^ l{Ei + E2), follows Wi f 1 and W2 -fa. 

For illustration, the cross section calculated with two neighbouring resonance states is 
shown in Fig. 1. The width of one of the states is fixed to Fi = 0.05 (in arbitrary units) 
while that of the other one is varied between F2 = 0.01 and 5.0. When F2 ^ Fi, it is F2 ~ 
const in the energy region of the narrow resonance, and the cross section shows a dip at the 
energy Ei. In this case, the broad resonance plays the role of a background for the narrow 
resonance. When the widths of both states are equal, the structure of the cross section is 
similar to that caused by (21) at the double pole of the S matrix. When F2 < Fi, the two 
peaks in the cross section are no longer symmetrical in relation to the center a.t E = 8. 

The corresponding couphng coefficients Wi and W2 (Figs. 2 and 3) show a resonance-hke 
behaviour at = {Ei + E2)/2. While the absolute values show the same tendency for both 
states, their phases behave differently. The phase of the broad state is almost not influenced 
by the interaction with the narrow one, while the phase of the narrow state jumps by 27r 
at the energy E = {Ei + E2)/2. This is caused by the minimum of Re(iyfc) that is reached 
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at E — (Ek + Ei)/2. The minimum value is — F; ?a F^ for the broad state, but — F; for 
the narrow one. In both cases, lm(Wk(i)) oscillates and vanishes at £^ = (E^. + Ei)/2. These 
results are in agreement with the fact that the broad resonance state plays the role of a 
background for the narrow resonance state. 

Of special interest is the case Fi = The resonance behaviour of the coupling coeffi- 
cients appearing at = (£"1 + E2)/2, is independently of the distance \Ei — E2\ of the two 
resonance states. It reflects the properties of a double pole of the S matrix. One clearly sees 
the phase jump by tt at = {Ei + E2)/2 (Figs. 2 and 3). This energy, E'^^ = {Ei + i?2)/2, 
is the critical value at which the wave functions of the two states are exchanged when the 
energy is parametrically varied (see Section IV). Here, W — Q. However, the phase jump of 
TT does not appear at every zero of the cross section, see Figs. 1 to 3. It follows further that 
the resonance behaviour of the coupling coefficients plays a role only for resonance states 
lying near to one another. 

As can be seen in Fig. 1, the interferences between the different resonance states cause, 
in the one-channel case, a separation of the peaks in the cross section. In this manner, the 
interferences may feign the existence of well isolated resonance states in spite of their strong 
overlapping. An extreme case are the two separated peaks appearing in the cross section 
with a double pole of the S matrix [2,10,14] or in its neighbourhood. Fig. 1. 

The line shape of the peaks in the cross section is described usually by means of energy 
independent Fano parameters. A recent example are the experimentally observed narrow 
peaks in the conductance through a quantum dot controlled by varying the strength of the 
magnetic field [22]. The energy independent Fano parameters are related to a representation 
of the S matrix (19) with energy independent [instead of the energy dependent in 
(20)]: 

Wl^tJl-i-~ ~ ,~ (22) 

^ \ Ek-Ei- |(F,-FO J ^ ' 

In the literature, mostly the W^. are used since they provide the standard parametrization 
of the S matrix. A recent example is the analysis of the bi-orthogonality of resonance wave 
functions in the molecule NO2 [23]. The representation (22) is equivalent to (20), except 
in approaching a double pole of the S matrix where (22) has a singularity in contrast to 
(20). The Wl lose, however, their physical meaning as coupling matrix elements between 
resonance states and continuum in the overlapping regime [10] . For this reason, we will not 
consider them in this paper. Instead, we will use the W}. that are meaningful also around 
double poles of the S matrix according to (21). 

It is easy to generalize the study to more than two resonance states. Suppose 

^ ^ TT fL_^ iili (23) 

^\ E — E + -F 

n=l ^ ~ 2 " 

instead of (17). In analogy with (19) and (20), it follows 

5 ^ 1 _ i V J^" . - (24) 

^ E — E ^ - V 

with 
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Y ( ^ ^ r'rn ' \ (25) 

\ Xn + Xjn XnXjji + XnXi + X^Xi J 

and Xn = E — En + fFn- The sum in the second term of (25) is running over m ^ n and 
that in the third term over m ^ n and / 7^ m, n. The denominator X„ + X^ of the second 
term is hnear in E while that of the third term, XnXm + XnXi + X^Xi, is quadratic in 
In any case, the Wn depend on energy in a non-trivial manner. 

As follows from (24) and (25), the couphng coefficients at a triple pole of the 5" matrix 
are Wn — r„/3 (in contrast to Wn = at a double pole). Here, the S matrix contains terms 
up to third order. When one of the widths is much larger than the other ones in the three- 
resonance case, Ti ^ (A; = 11,171), it is Wi ~ Ti and Wk = — when E = ^{Ei + E^). 
These relations are in complete analogy with those obtained for the two-resonance case. 

In Figs. 4 and 5, the cross sections with three resonance states are shown, two of which 
are lying symmetrically around the position of the third one at = £'3 = 8, as well as the 
coupling coefficients W-^. The different curves in Fig. 4 are obtained by varying the width 
F3 between 0.05 and 5. In the first case, the widths of all three states are equal while in 
the other cases, the middle state overlaps the narrower ones. When F3 ^ F^, {k = 1,2), 
the broad state can be considered as a "background" for the two other ones: they appear 
as dips in the cross section. This can be seen better in Fig. 6 where one broad resonance 
state (r3 = 3.0) is shown that overlaps two narrow ones (Fi = ^2 = 0.05). In both cases, 
a peak appears in the middle of the spectrum at = 8 in contrast to the two-resonance 
case Fig. 1. In Fig. 5, the widths of all three states are kept constant while the distance 
between them is varied. Altough the peaks in the cross section seem to be well isolated from 
one another for (Ek — Ek±i)/(Tk -\- Ffeii) = 2.5 (full curve), the interferences between the 
resonance states are not vanishing. 

Interesting are again the phase jumps of the coupling coefficients 1^3 at the critical values 
of the energy E'^'^ = [E^ + E^)/2; k = 1,2. They are mostly smaller than those of Wi and 
W2 (not shown). The influence of a broad resonance state on the coupling coefficients of 
narrow ones is illustrated in Figs. 6 and 7 with the state 3 being, respectively, much broader 
and much narrower than the states 1 and 2. At the energies 6 and 9, the states 1 and 2, 
respectively, interact with the state 3 while the interference picture at the energy 7 is caused 
by the two states 1 and 2 with equal widths. The phase jump of tt at this energy (Fig. 7) 
is reduced in the presence of the broad resonance state (Fig. 6). Also ll^il and 11^21 differ 
in the two cases without and with a broad state. 

Measurements of phase and magnitude of the reflection and transmission coefficients of 
a quantum dot are performed in [24]. As a result, the phases of the dot's transmission and 
reflection coefficients change abruptly by about tt at some energy in the resonance peak. 
The phase changes are very similar to those observed in the present calculations for the 
W in the one-channel case (see the figures). As discussed above, they are related to the 
unitarity of the S matrix in the overlapping regime. The results are expected to be similar 
for the case with more channels (or terminals) since, as will be shown in Sect. IV, they are 
characteristic of the intrinsic wave functions and ^f^, respectively, of the system. 

In the consideration presented in the present section, the E^ and F^ are assumed to be 
independent of the energy. This is, mostly, a good approximation in the energy range of 
the resonance state k. In any case, the energy dependence of the coupling coefficients Wk 
arises primarily from their resonance behaviour caused by a neighbouring resonance state. 
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Eq. (20). It may be influenced, of course, by the energy dependence of and F^, especially 
when the levels repel or attract each other by varying a certain parameter. 

IV. WAVE FUNCTIONS NEAR AVOIDED LEVEL CROSSINGS IN THE 

COMPLEX PLANE 

The couphng coefficients between system and continuum are defined by (7^)^, Eq. (5). 
Their energy dependence and phase are determined, in the one-channel case, by the energy 
dependence and phase of the (after removing the common phase and energy depen- 

dence caused by the C,e)- Much can, therefore, be learned on the behaviour of the coupling 
coefficients between system and environment from a study of the wave functions ^k- 

Characteristic of the overlapping regime are avoided level crossings in the complex plane. 
At an avoided level crossing, the wave functions of the two crossing states are exchanged. 
This fact is very well known for a long time for discrete states (Landau-Zener effect). It 
holds also for resonance states in the adiabatic limit [2]. The difference between avoided 
crossings of discrete and resonance states consists mainly in the fact that resonance states 
may cross in the complex plane even when the interaction between them is non-vanishing 
and, furthermore, that the crossing may lead not only to level repulsion but also to level 
attraction. While level repulsion is accompanied by the tendency to equilibrate the widths 
of the resonance states, level attraction is accompanied by the formation of different time 
scales in the system. 

The avoided level crossing of resonance states can be traced back to a branch point in 
the complex plane [25]. When the conditions for crossing of the two levels are fulfilled, the 
branch point in the complex plane is nothing else than a double pole of the S matrix (see 
also [26]). In any case, the wave functions of the two levels are exchanged at the critical 
value of the parameter at which the two levels avoid crossing (or cross in one point). Here, 
either the widths or the energies (positions) of both states are equal. In the first case, the 
avoided crossing happens, as for discrete states, in the energies of the resonance states traced 
as a function of the considered parameter. In the second case, however, it appears in their 
widths. 

Further studies have shown that the wave functions at the critical value a"^^ of the pa- 
rameter a are exchanged according to 

h^±i^i (26) 

where k and I are the two crossing states. This relation has been obtained analytically [25] 
as well as in a numerical study of atoms in a laser field [16]. It is related to nonlinear terms 
that appear in the Schrodinger equation due to the biorthogonality of the eigenfunctions of 
the effective Hamiltonian Ti.. At a'^^, the sign of the imaginary part of the wave function $fe 
jumps from -|- to — (or opposite) even when the two states avoid crossing and = \^k\ 
remains finite [25]. That means, in a certain parameter range a™™ < a^'^ < a™^, the wave 
functions of the two states k and I ^ k are mixed, 

^f = f3kh±i(3i^i. (27) 
The wave functions change smoothly (without any jump of the sign of their components) 



12 



from (3k 
to (3k 



±1, (3i^0 at a ^ a"^'" < a' 
0, A ^ ±1 at a ^ a""^ > a' 



cr 



(28) 



The values a™™ and a™^ may be quite different from one another [25]. Only in the case the 
avoided level crossing shrinks to one point, being the double pole of the 5" matrix, (3k — or 
±1 for all a but a'^^. In any case, the two wave functions {^k ± i ^i) and -^{^i^i ^k) 
remain unchanged at a = a" under the exchange (26). 

Hence, in the parameter range a""° < a < a™^^, the wave functions of the two states are 
^f", but not $fc {k = 1, 2). The two wave functions are restored, after the exchange at a", 
only at a > a™^. In other words, using the representation 



9k depends on a when a™™ < a < a^^. Beyond this parameter area, 9k depends on a much 
weaker (if at all). After removing a common phase factor, it follows from (28) for an avoided 
level crossing near the double pole: 6k ±7r/4 and ±37r/4, respectively, in approaching a'^'^ 
and ~^ or TT in approaching a™'" or a^^^. At a'^'" the states are chiral. 

The comparison of (26) with experimental data obtained from microwave cavities has 
been discussed in detail in [27]. All the data published in [28,29] for the case of two levels 
of the system that are well isolated from the other ones, can be explained by means of (26). 
Some chirality appears: left and right turns around the double pole of the S matrix with 
(3k = (3i are different from one another according to (27) and (28). 

Phase changes of the wave functions in the conductance through a microwave cavity 
have been considered in [30,31]. The question is studied to what extent the transport 
through the cavity changes the structure of its internal wave functions. It is demonstrated 
theoretically as well as experimentally [30] that the standing waves of the original cavity 
are transformed more or less completely into running waves propagating from the entrance 
antenna to different exit ports. This is expressed by the two hmiting cases: (i) the real and 
imaginary parts of the wave functions are strongly correlated, and (ii) they are completely 
uncorrelatcd. In the last case, the real and imaginary parts of the wave functions of the 
resonance states evolve independently in the open microwave cavity. 

In [32], the phase difference between two modes has been measured in a cavity composed 
of two almost identical semi-circular parts. The two modes are each localized in one of 
the semi-circular parts of the cavity and are excited separately by appropriately positioned 
dipole antennas. The corresponding two eigenvalues are well separated from all the other 
ones. By varying two parameters (designed here together by a) of the cavity, their avoided 
crossing in the complex plane to an (almost) true crossing can be traced. Furthermore, the 
eigenfunctions are studied by mapping the distributions of the electric field. Finally, the 
phase difference A — 9k — 9i between the antennas has been found for different distances of 
a from the critical value a". The results obtained in the experiment [32] are A ^ 7r/2 for 
a a'^'^ and A = tt for a beyond the range a™™ < a < a^'^^. They agree with (27) and (28) 
by using the representation (29) [33]. 

According to (5), the phase of the coupling coefficients (7^:)^ is determined by that of 
i^kY- Considered as a function of a certain parameter, the phases of both expressions vary, 
in the one-channel case, in the same manner. Since the parameter may be also the energy of 
the system [2], the phase of the (7^)^ varies generically with the energy in the same manner 




(29) 
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as the phase of the ($fc)^. The results discussed above are in full agreement with those 
following immediately from the unitarity of the S matrix, Figs. 2 to 7. Even the jumps by tt 
appearing in the phases of the coupling coefficients W in Figs. 2, 3 and 7 can be explained 
by Eqs. (26) and (27). In other words, Eqs. (26) to (28) coincide with the postulation 
of the unitarity of the S matrix. In any case, nonlinearities are responsible for the energy 
dependence of the coupling coefficients between system and continuum. 

In realistic systems, mostly more than two resonance states are coming close to one 
another, i.e. avoided crossings of more than two resonance states take place in the parameter 
range a™™ < a < a™*^^. As an extreme case, the branch point in the complex plane may 
be a multiple pole of the S matrix. The relation between the different wave functions is, 
in such a case, more complicated than that for two states where the wave functions of only 
two resonance states are exchanged according to (26), mixed in the range a™" < a < 
and restored beyond this range. Such a situation is studied experimentally [28] as well 
as theoretically in different approaches [4,25,34,35], compare also Figs. 2 and 3 for two 
resonance states with Figs. 4 and 5 for three resonance states as well as Fig. 6 with Fig. 7. 

V. CONCLUDING REMARKS 

The resonance phenomena are described well by two ingredients also at high level density. 
The first ingredient is the efi^ective Hamiltonian Ti. that contains all the basic structure 
information involved in the Hamiltonian H, i.e. in the Hamiltonian of the corresponding 
closed system with discrete eigenstates. Moreover, Ti, contains the coupling matrix elements 
W^' between discrete states and continuum that account for the changes of the system 
under the influence of its coupling to the continuum. These matrix elements arc responsible 
for the non-Hermiticity of Ti. and its complex eigenvalues which transfer the discrete states 
into resonance states and determine not only their positions but also their (finite) lifetimes. 

The second ingredient is the unitarity of the S matrix that has to be fulfilled in all 
calculations of resonance phenomena. The unitarity of the S matrix causes a non-trivial 
energy dependence of the coupling matrix elements between resonance states and con- 
tinuum. This energy dependence becomes decisive in the overlapping regime even in the 
case the lifetimes of the overlapping states are very different from one another and the dif- 
ferent long-lived states seem to be well isolated from one another. It is taken into account 
in the unified description of structure and reaction aspects since any statistical or pertur- 
bative assumptions are avoided in solving the basic equation (1). The unitarity of the S 
matrix influences also the phases of the wave functions of the resonance states that change 
generically in approaching avoided level crossings in the complex plane. 

In the non-overlapping regime, both ingredients are fulfilled in almost all theoretical 
approaches. Here, the wave functions and positions of the resonance states are described, 
to a good approximation, by the wave functions and positions of the discrete states of the 
corresponding closed system. The coupling matrix elements between system and continuum 
can be calculated by means of the wave functions of the discrete states. They are energy 
independent, to a good approximation. 

In the overlapping regime however with many avoided level crossings, the wave functions 
of the resonance states suffer phase changes. Furthermore, the coupling coefficients of the 
resonance states to the continuum show a resonance-like behaviour caused by the interaction 
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with a neighboured resonance state. Both effects are described by nonhnear terms appearing 
in the Schrodinger equation and the S matrix, respectively. They are related to one another 
and can not be neglected at high level density. 

As a conclusion of these results, the coupling of a quantum system to the environment 
may change its properties. The changes are small as long as the coupling strength between 
system and environment is smaller than the distance between the individual states of the un- 
perturbed system, i.e. smaller than the distance between the cigenstatcs of the Hamiltonian 
H. The changes can, however, not be neglected when the coupling to the continuum is of 
the same order of magnitude as the level distance or larger. In such a case, the changes can 
be described neither by perturbation theory nor by introducing statistical assumptions for 
the level distribution. Here, nonlinear effects become important which cause a redistribution 
of the spectroscopic properties of the system and, consequently, changes of its features. 

Under the influence of the coupling to the continuum, not only level repulsion but also 
level attraction may appear that are accompanied by the tendency to form a uniform time 
scale for the system in the first case, but different time scales in the second case. The 
formation of different time scales in an open quantum system that is accompanied by level 
attraction, is accompanied also by the appearance of a non-trivial energy and phase de- 
pendence of the coupling coefficients W. The use of an effective non-Hermitian Hamilton 
operator in describing scattering processes in the overlapping regime is therefore meaningful 
only when, at the same time, the energy dependence of the W is considered. 

Acknowledgment: 1 am indebted to A.I. Magunov, N. Moiseyev and E. Persson for valu- 
able discussions and to H. Schomerus for a critical reading of the manuscript. 



15 



REFERENCES 



[1] H. Feshbach, Ann. Phys. (N.Y.) 5, 357 (1958) and 19, 287 (1962). 

[2] J. Okolowicz, M. Ploszajczak, and I. Rotter, Phys. Rep. 374, 271 (2003). 

[3] Y.V. Fyodorov and B. Mehlig, Phys. Rev. E 66, 045202 (2002); Y.V. Fyodorov and 

H.J. Sommers, J. Phys. A 36, 3303 (2003). 
[4] C. Jung, M. Miiller and I. Rotter, Phys. Rev. E 60, 114 (1999). 
[5] T. Guhr, A. Miiller- Groehng, and H. A. Weidenmiiller, Phys. Rep. 299, 190 (1998). 
[6] K. Pichugin, H. Schanz, and P. Seba, Phys. Rev. E 64, 056227 (2001). 
[7] Y.V. Fyodorov and H.J. Sommers, J. Math. Phys. 38, 1918 (1997). 
[8] P.M. Dittes, Phys. Reports bf 339, 215 (2000). 

[9] D.V. Savin, V.V. Sokolov, and H.J. Sommers, Phys. Rev. E 67, 026215 (2003). 
[10] A.I. Magunov, I. Rotter and S.I. Strakhova, to be published 
[11] I. Rotter, Rep. Prog. Phys 54, 635 (1991). 

[12] E. Persson, I. Rotter, H.J. Stockmann, and M. Earth, Phys. Rev. Lett. 85, 2478 (2000); 
H.J. Stockmann, E. Persson, Y.H. Kim, M. Barth, U. Kuhl, and I. Rotter, Phys. Rev. 
E 65, 066211 (2002). 

[13] R.G. Newton, Scattering Theory of Waves and Particles, Springer- Verlag New York, 
1982. 

[14] M. Miiller, P.M. Dittes, W. Iskra and I. Rotter, Phys. Rev. E 52, 5961 (1995). 
[15] C.J. Joachain, M. Dorr, and N. Kylstra, Adv. At. Mol. Opt. Phys. 42, 225 (2000). 
[16] A.I. Magunov, I. Rotter and S.I. Strakhova, J. Phys. B 32, 1669 (1999); J. Phys. B 34, 
29 (2001). 

[17] S.D. Bosanac, Phys. Rev. A 30, 142 (1984) and Phys. Rev. A 30, 148 (1984). 

[18] T. Taniguchi and M. Biittiker, Phys. Rev. B 60, 13814 (1999). 

[19] E. Hernandez, A. Jauregui and A. Mondragon, J. Phys. A 33, 4507 (2000). 

[20] W. Vanroose, Phys. Rev. A 64, 062708 (2001). 

[21] A. Silva, Y. Oreg and Y. Gefen, Phys. Rev. B 66, 195316 (2002). 

[22] V. Madhavan, W. Chen, T. Jamneala, M.F. Crommie, and N.S. Wingreen, Science 280, 
567 (1998); J. Gores, D. Goldhaber-Gordon, S. Heemeyer, M. A. Kastner, H. Shtrikman, 
D. Mahalu, and U. Meirav, Phys. Rev. B 62, 2188 (2000); K. Kobayashi, H. Aikawa, S. 
Katsumoto, and Y. lye, Phys. Rev. Lett. 88, 256806 (2002). 

[23] S.Yu. Grebenshchikov, J. Phys. Chem. A 104, 10409 (2000). 

[24] A. Yacoby, M. Heiblum, D. Mahalu, and H. Shtrikman, Phys. Rev. Lett. 74, 4047 
(1995); E. Buks, R. Schuster, M. Heiblum, D. Mahalu, V. Umansky, and H. Shtrikman, 
Phys. Rev. Lett. 77, 4664 (1996); R. Schuster, E. Buks, M. Heiblum, D. Mahalu, V. 
Umansky, and H. Shtrikman, Nature (London) 385, 417 (1997). 

[25] I. Rotter, Phys. Rev. E 64, 036213 (2001); Phys. Rev. C 64, 034301 (2001). 

[26] N. Moiseyev, Phys. Rep. 302, 211 (1998). 

[27] I. Rotter, Phys. Rev. E 65, 026217 (2002); Phys. Rev. E 67, 026204 (2003). 

[28] H.M. Laubcr, P. Wcidcnhammcr, and D. Dubbers, Phys. Rev. Lett. 72, 1004 (1994). 

[29] C. Dembowski, H.D. Graf, H.L. Harney, A. Heine, W.D. Heiss, H. Rehfeld, and A. 

Richter, Phys. Rev. Lett. 86, 787 (2001). 
[30] P. Seba, F. Haake, M. Kiis, M. Barth, U. Kuhl, and H.J. Stockmann, Phys. Rev. E 56, 

2680 (1997). 



16 



[31] H. Ishio, A.I. Saichev, A.F. Sadreev, and K.F. Berggren, Phys. Rev. E 64, 056208 
(2001). 

[32] C. Dembowski, B. Dietz, H.D. Graf, H.L. Harney, A. Heine, W.D. Heiss, and A. Richter, 
Phys. Rev. Lett. 90, 034101 (2003). 

[33] The interpretation of the experimental data given here on the basis of Eqs. (27) and 
(28) differs from that given in [32]. The authors of [32] consider the reduction of the 
dimension of the Hilbert space at a = a'^^ (" defect of the Hilbcrt space" ) that is caused 
by the coalescence of two eigenmodes. In the case studied, the number of states is 
assumed to be reduced from 2 to 1 at a^^ . This reduction of the number of states is 
suggested to be responsible for the reduction of A from tt to 7r/2 at a^^ . In difference 
to that interpretation, the reduction of A from tt to 7r/2 is caused, in the formalism 
given in the present paper, first of all by the different degree of mixing of the two states 
according to (28). It takes place smoothly in a certain region around a™ in agreement 
with the experimental data published in [32] . In the present formalism, diabolic points 
and exceptional points are related to one another. Nevertheless, they are not the same 
in contrast to the statement given in reference [5] of [32]. 

[34] F. Pistolesi and N. Manini, Phys. Rev. Lett. 85, 1585 (2000); N. Manini and F. Pistolesi, 
Phys. Rev. Lett. 85, 3067 (2000). 

[35] E. Narevicius and N. Moiseyev, Phys. Rev. Lett. 84, 1681 (2000); J. Chem. Phys. 113, 
6088 (2000). 



17 



FIGURES 




FIG. 1. Cross section with two resonance states at Ei = 7.99 and E2 = 8.01. The width of 
one of the states is fixed to f 1 = 0.05, while that of the other state is varied: f 2 = 5.0 (dashed 
curve), r2 = 0.05 (full curve), r2 = 0.01 (dash-dotted curve). The dotted curves are calculated 
with r2 = 1-0, 0.5, 0.1, 0.025, respectively. Cross section and energy are given in arbitrary units. 



18 




19 




41 < < < 1 

7.8 7.9 8 8.1 8.2 

E 



FIG. 3. Coupling coefficient W2 of the resonance state at E2 = 8.01. Tlie position and widtfi 
of tfie otfier state as well as the width r2 is as in Fig. 1. In the upper part of the figure, \W2\ is 
multiplied by 1 (dashed curve), 100 (full curve) and 500 (dash-dotted curve) and by 5, 10, 50, and 
200, respectively (dotted curves). 
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FIG. 4. Cross section (top) with three resonance states at Ei = 7.99, E2 = 8.01, = 8.0 and 
fi = fs = 0.05, fa = 0.05 (full curve), fg = 1.0 (dash-dotted curve), and fg = 5.0 (dashed curve). 

Coupling coefficient (middle and bottom) IV3 of the resonance state in the middle of the spectrum. 
IW3I is multiphed by 100 (full curve), 5 (dash-dotted curve) and 1 (dashed curve). Cross section 
and energy are given in arbitrary units. 



21 




FIG. 5. Cross section (top) with three resonance states of widths Fi = r2 = F3 = 0.05 at 
El = 7.99, E2 = 8.01, £"3 = 8.0 (dotted curve), Ei = 7.9, E2 = 8.1, E^ = 8.0 (dash-dotted 
curve). El = 7.75, E2 = 8.25, E^ = 8.0 (full curve). Couphng coefficient W3 (middle and bottom) 
of the resonance state in the middle of the spectrum. Cross section and energy are given in arbitrary 
units. 
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FIG. 6. Cross section (top) with one resonance state of width T3 = 3.0 at = 8.0 (dotted 
curve) and with two additional resonance states of widths Fi = r2 = 0.05 at Ei = 4.0 and 
E2 = 10.0, respectively (full curve). Coupling coefficients (middle and bottom) Wi (dashed), W2 
(dash-dotted) and W3 (full curve) of the three resonance states. \Wi\ and \W2\ are multiplied by 
60. Cross section and energy are given in arbitrary units. Note the different energy scales for the 
cross section and the coupling coefficients. 
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FIG. 7. The same as Fig. 6, but fg = 3 • 10"^ 
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